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Abstract: In view of precision studies of the Higgs sector at the Run II of the LHC, the 
improvement of the accuracy of the theoretical prediction is becoming a pressing issue. In 
this framework, we detail a calculation of the full Next-to-Leading Order (NLO) electroweak 
corrections to Higgs boson decay into four charged leptons, by considering the gold-plated 
channel H —)■ Z^*'^ —)• 2121', l,l! = e,[i. We match the NLO corrections with a QED 
Parton Shower (PS), in order to simulate exclusive multiple photon emission and provide 
novel results at NLOPS electroweak accuracy. We compare our NLO predictions to those of 
the program Prophecydf and present NLOPS phenomenological results relevant for Higgs 
physics studies, with particular attention to precision measurements of the Higgs boson 
mass, spin-parity assignment and tests of the Standard Model. Our calculation is imple¬ 
mented in a new code, Hto41, which can be easily interfaced to any generator describing 
Higgs boson production. As an example, we provide illustrative results for Higgs produc¬ 
tion and decay in the process gg ^ H ^ M using POWHEG with NLOPS accuracy in the 
production mode. 
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1 Introduction 

With the announcement in 2012 of the discovery of a new particle in the search for the 
Standard Model (SM) Higgs boson by the ATLAS [1] and CMS [2] collaborations at the 
CERN LHC, particle physics entered a new era. The data collected at the centre-of-mass 
(c.m.) energies of 7 and 8 TeV have been analyzed by the two experiments in order to 
establish whether the newly discovered particle is actually the boson predicted in the SM 
as relic of the mechanism of electroweak symmetry breaking (EWSB) [3-8]. 

The mass of the observed particle has been precisely measured by studying the two 
cleanest decay channels given by the decays into a photon pair and into four charged 
leptons. The combination of the two channels —?■ 77 and 77 —?■ (4e, 4/i, 2e2/i), which 
have excellent mass resolution and where excesses with large significance are observed [9- 
14], presently provides a mass measurement of approximately 125 GeV for each experiment, 
with a relative uncertainty of better than 0 . 2 % for the combined ATLAS-CMS measurement. 

Concerning the main production mechanisms of the SM Higgs boson at hadron colliders, 
i.e. gluon-gluon fusion, vector boson fusion (VBF), associated production with a massive 
vector boson and associated production with top quarks, the studies performed at the LHC, 
based on the analysis of individual production signal strengths for various decay modes, have 
provided a clear observation of Higgs production through gluon fusion and an evidence for 
VBF production, with a significance above the 3a level, and for associated VH{V = W,Z) 
production at about 3a [12, 15]. 

Various tests of the couplings of the new particle to bosons and fermions have been 
carried out both by ATLAS and CMS collaborations. In particular, the measured ratio 
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of the couplings of the Higgs particle to W and Z bosons, which is an important probe 
of the EWSB mechanism as fixed by the custodial symmetry, is compatible with the SM 
expectation and, more generally, no significant deviation from the SM is observed from the 
coupling strength studies [9, 12, 15]. Noticeably, evidence for the direct coupling of the 
Higgs boson to down-type fermions has been reported through the study of the challenging 
decay modes Higgs into bottom quarks and r leptons [16, 17]. 

Last but not least, the spin and parity quantum numbers of the discovered particle 
have been assessed, by means of a systematic analysis of its production and decay processes. 
The data strongly favor the scalar nature = 0“’' of the observed particle, while rejecting 
other non-standard hypotheses (J^ = 0“,1^,2'’') or possibility of CP mixtures at high 
confidence [14, 18, 19]. 

All these measurements marked the start of a new era of precision Higgs physics and 
were accompanied by an impressive theoretical effort summarized in three CERN reports 
by the LHC Higgs Cross Section Working Group [20-22]. These studies, as well as the 
related theoretical work, are in continuous progress and will continue during the Run H of 
the LHC at higher energies and luminosity. 

In this paper, we focus on the Higgs boson decay into four charged leptons, i.e. 77 —t- 
in order to provide novel precision predictions of interest for future studies 
of the Higgs sector at the LHC. This decay channel plays a particularly relevant role, as 
it provides the cleanest experimental signature, given by a peak in the four lepton mass 
spectrum on top of a flat and small background. Actually, the 77 —)• 47 decay mode allows 
to derive a precise mass measurement in the different combinations of lepton final states, 
to assess the spin-parity quantum numbers using sensitive angular distributions and to 
perform precision tests of the SM at the level of differential cross sections [23]. In the 
off-shell region, the 77 —>■ 47 data can also be used to put constraints on the total width of 
the Higgs boson [24, 25]. 

In the light of the above motivations, we compute the full set of next-to-leading or¬ 
der (NLO) electroweak corrections to 77 —)■ Z^*'> —)• 47, with 47 = 4e, 4/r, 2e2/r. We 

match the NLO corrections to a QED Parton Shower (PS), in order to simulate multi¬ 
ple photon emission exclusively and provide final results at NLOPS electroweak accuracy. 
The calculation is available in an event generator, Hto41^, which can be interfaced to any 
code describing Higgs boson production. The PS approach is based on the ideas first 
presented in Ref. [26] for the simulation of the Bhabha scattering process at GeV-scale 
e'’'e“ colliders and later applied to the study of single WjZ production in hadronic col¬ 
lisions [27, 28]. The matching procedure is a generalization of the method developed in 
Refs. [29, 30] for the precision calculation of 2 —2 processes in QED (as encoded in 
the program BabaYagaONLO [31, 32]) and also implemented in the event generator Horace 
for the calculation of single WjZ hadroproduction processes at NLOPS electroweak accu¬ 
racy [33, 34]. 

The NLO electroweak and QCD corrections to 77 —t- 4 fermions decay processes have 
been calculated in Refs. [35, 36] and are available in the Monte Carlo (MC) program 

^The reference web page is http:/www.pv.infn.it/hepcomplex/hto41 .html 
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Prophecy4f [37, 38], which is used in the context of Higgs studies at the LHC for the 
precision calculation of the branching ratios of the decays H —?■ —)• 

4 fermions. In Prophecy4f higher-order photonic corrections are taken into account in 
terms of QED collinear Structure Functions. A preliminary study of the impact of the 
gauge-invariant NLO QED and PS corrections to the determination of the Higgs boson 
mass in the H ^ M decay was performed in Refs. [39, 40]. 

The article is organized as follows. In Section 2 we describe the details of our calcula¬ 
tion, with particular emphasis on the method used for the matching of the NLO electroweak 
corrections with the QED PS. In Section 3 we present our phenomenological results: in Sec¬ 
tion 3.1 we show a sample of comparisons between our predictions and those of Prophecy4f 
as a benchmark of the NLO computation, in Section 3.2 we provide results for various 
observables at NLOPS EW accuracy, while in Section 3.3 we present the results for Higgs 
production and decay in the channel gg ^ H ^ Ai obtained in terms of POWHEG [41, 42] 
interfaced to Hto41. In Section 4 we draw our conclusions. 

2 Details of the calculation 

2.1 Next-to-leading order (NLO) electroweak corrections 

The NLO electroweak corrections to the Higgs boson decay into four charged leptons consist 
of QED and purely weak contributions. Since the 77 —)• 4£ decay is a neutral-current process, 
the two subsets are separately gauge invariant and can be computed separately as well. 

The 0{a) QED corrections are obtained by attaching a virtual or real photon to each 
charged lepton leg. They are expected a priori to provide the dominant contribution, as 
photons which are emitted collinear to a lepton give rise to large logarithmic corrections of 
the form oTog (m^/Q^), where me is the lepton mass and Q some typical energy scale. 

The QED virtual corrections comprise vertex and pentagon diagrams (in the on-shell 
renormalization scheme), while real photon corrections are induced by the bremsstrahlung 
process H —>- 47 - 1 - 7 . The two contributions are separately infrared (IR) divergent but their 
sum is IR-finite. We treat the IR singularity according to the standard QED procedure 
of assigning a small fictitious mass A to the photon in the computation of the virtual and 
real contributions. More precisely, the Higgs decay width associated to the bremsstrahlung 
correction is separated in two pieces and calculated as follows (in a shorthand notation) 

J dTreai = J d4>5 |A4o(i7 -> 47) X (eikonal factor) 

+ d<I>5|W((77^47 + 7)P ( 2 . 1 ) 

where Mh is the Higgs mass, e is a soft-hard energy separator (e <C Mh), —> 47) is 

the amplitude of the lowest-order (LO) process 77 —> 47 and Af (77 —> 47 -|- 7 ) is the matrix 
element of the radiative decay process 77 —> 47 -|- 7 , d 4>5 being the 4 leptons plus 1 photon 
phase space element. In Eq. (2.1) eikonal factor stands for the analytical expression of the 
real radiation correction in the soft limit —> 0. The integral in the first line can be done 
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analytically (see e.g. [43]) and the one in the second line is performed using standard MC 
integration with importance sampling. 

The QED virtual counterpart is computed according to the following formula 

dr^im = ^d4>42Re{Af?'^°(A1o)*} (2.2) 

where is the one-loop amplitude associated to the 0{a) vertex and pentagon dia¬ 

grams. 

We perform the IR cancellation by taking the sum of Eq. (2.1) and Eq. (2.2) in the 
numerical limit A —?■ 0. As a cross-check of the calculation, we tested that the inclusive NLO 
QED correction coincides with 2-3/4 (a/vr), which is correctly twice the inclusive final-state 
0(a) electromagnetic correction to the Z —?• decay [44]. An important comment is 

in order here. The tree-level amplitude, as well as the amplitude for the real radiation 
process, contains poles in the phase space, corresponding to the points where the momenta 
of the £~^£~ pairs and of the (£~^£~'y) system cross the zero of the inverse propagators: 
{P£+ +P£-)^ = or {P£+ + P£- +4*7)^ = ^z- These poles are avoided considering that 
the Z boson is an unstable particle, i.e. its propagator contains the finite Z-width. This, 
however, would spoil the IR cancellation between real and virtual corrections of Eq. (2.1) 
and Eq. (2.2), respectively, unless in Eq. (2.2) the QED virtual corrections are calculated 
with unstable Z bosons. The scheme which we adopt for the introduction of the width in 
the Z boson propagator, without introducing gauge invariance violations, is the complex 
mass scheme [45, 46], which also allows us to include weak loop corrections consistently.^ 

Concerning the basic features underlying the computation of the complete 0(a) virtual 
corrections, we briefly describe the most important aspects in the following. Since we work 
in the’t Hooft-Feynman gauge, all the particles present in the spectrum of the SM, including 
the Fadeev-Popov and Higgs-Kibble ghosts, are involved in the calculation. The correspond¬ 
ing Feynman diagrams include, in addition to two-point functions, rank-two tensor three-, 
four- and five-point functions. The related ultraviolet divergencies are regularized by means 
of dimensional regularization. The reduction of the tensor n-point functions is carried out 
by means of the symbolic manipulation program FORM [47]. The necessary scalar form fac¬ 
tors with complex masses are evaluated using Looptoolsv2.10 [48, 49], which implements 
the evaluation of the reduction of tensor five-point integrals according to Refs. [50, 51], 
as well as according to Passarino-Veltman reduction techniques [52]. The form factors are 
calculated with complex masses and real external squared momenta. This is sufficient for 
the implementation of the “simplified version of the complex renormalization”, as described 
in Refs. [45, 46]. The complete expressions for the counterterms in the on-shell scheme and 
for the basic self-energy diagrams are taken from Ref. [53]. Since the collinear singularities 
associated to the photon becoming collinear with one of the leptons are regulated by the 
finite lepton mass, the kinematics of the radiative process is calculated including exactly 
the contribution of lepton masses. In order to allow the cancellation of soft IR singularities, 

^Actually, the complex mass scheme is used in our calculation of the weak contributions due to the 
exchange of W bosons and top quarks as well, where a complex top mass is introduced, in particular, to 
evaluate loop diagrams with internal top quarks when the Higgs mass is close to the tt threshold. 
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also the tree-level kinematics is calculated with complete lepton mass effects taken into ac¬ 
count. In addition, this gives automatically the correct phase space integration boundaries 
for the diagrams of the virtual contribution where a virtual photon is connected to one 
external lepton pair. Although the kinematics is treated exactly, the non-IR 0{a) virtual 
amplitudes are calculated in the approximation of neglecting finite fermion mass effects 
(with the exception of the quark Yukawa couplings, e.g. in the fermion-loop Higgs vertex 
corrections). These contributions are neglected in our calculation as they are irrelevant in 
view of a target theoretical accuracy of the order of 0 . 1 % and their inclusion would make 
the numerical computation more time consuming. 

In formulae, the Higgs width including one-loop weak corrections is obtained as 

= ^d^^2Re{Mr^HMor} (2.3) 

where is the one-loop amplitude associated to the full set of 0{a) weak diagrams. 

To check some relevant ingredients of our calculation of one-loop weak corrections, we 
compared our predictions for the Higgs decays i/ —)• ZZ, 77 at NLO electroweak accuracy 
with those of Ref. [54], finding perfect agreement. 

In conclusion, our predictions for the Higgs boson decay into four leptons at NLO EW 
accuracy are given by the sum of Eq. (2.1), Eq. (2.2) and Eq. (2.3), supplemented with the 
necessary renormalization conditions. 


2.2 Matching NLO electroweak corrections to QED Parton Shower 

In the present section, we sketch our scheme for the matching of the NLO EW corrections 
with a QED PS. We closely follow the approach already presented and successfully applied 
to QED processes at low energies and Drell-Yan WjZ production at hadron colliders [29, 
30, 33, 34]. 

On general grounds, the partial decay width corrected for the emission of an arbitrary 
number of photons in a PS framework can be written as follows: 

1 °° 1 

dr™ = ^n({p),e)^- |A<P({p).W)r d<f„({p},{t)) (2.4) 

n=0 

where {p,k} stands for the set of the final state lepton and photon momenta pi, p 2 , Ps, 
P 4 , fci, • • • , kn, (of order a”) is the PS approximation to the squared amplitude for 

the decay H ^ 4i + ny, d<I>„ is the exact phase space for the decay and n({p},e) is the 
Sudakov form factor accounting for unresolved emission, i.e. soft (up to a cut-off energy e) 
and virtual corrections in the PS approximation. It is understood that the integral over the 
phase space has a lower limit for the photon energies set to e, to ensure the cancellation of 
the IR divergencies. 

The quantities d4>„ and n({p},e) read explicitly 


d<I'n = 


(2vr)3n+8 


^(4) 


Ph - 


i=l 


TT TT 

li 9 ti0 11 2k^ 


2p' 

j=i i=i 


(2.5) 
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( 2 . 6 ) 


n(M,e) = exp 


-—Lie 

27r 


0\2 


L = / dJ2fc(A:") 


ij=l 


Pi ■ Pj 


{Pi ■ k)iPj ■ k) 


In Eq. (2.6), L generates the soft/virtual collinear logarithms, including also interferences 
effects of radiation coming from different charged legs, and le, the integral of the Altarelli- 
Parisi vertex for the branching £ —)> £ + 7 , generates the infrared logarithms. It is explicitly 
given by: 


Ie= [ =-21ne-^+2e- (2.7) 

Jo I-z 2 2 ^ ^ 

In the dehnition of L, the integral is performed over the angular variables of k, and rji equals 
1 if i is an anti-fermion or —1 if it is a fermion. 

The integral over the phase space as in Eq. (2.5) is performed after choosing a convenient 
set of independent variables and using multi-channel MC importance sampling techniques 
to improve the integration convergence and follow the peaking structure of the partial decay 
width of Eq. (2.4) to help events generation. The fully exclusive information on hnal state 
particles momenta is kept. Details of the implementation are given in appendix A. 

Before discussing the inclusion of NLO corrections into Eq. (2.4), it is interesting to 
point out that the squared amplitudes with photon emissions are enhanced in regions of 
the phase space where the photons are soft and/or collinear or where the Z propagators 
are resonating. In this perspective, a good approximation to the exact matrix element can 
be written in the form^: 

__A nPiPVi-^niki) ^2 8) 

^ (pi +P 2 + Qvf - (7*3 +Pi + Rvf - r=i 

In the previous equation, c is a shorthand for the HZZ coupling, {(T,t} label fermion 
and photon elicities, = Uai{Pi)Y{gv - gAlb)vaj{Pj), "i| = M| - iTzMz and £r{k) are 
the photon polarization vectors. "P is a n-dimensional vector whose component is the 
index of the fermion to which the photon is attached and the sum over V denotes all 
possible ways to share n photons among the four fermions. Finally, Q-p is the sum of the 
momenta of the photons, for a given V, attached to the electron current (R-p to the muon 
current). 

Equation (2.8) is derived from the amplitude for the emission of photons in the soft 
limit but keeping the dependence on the photon momenta in the Z propagators. The sum 
over the elicities of the squared amplitudes of Eq. (2.8) gives an approximation of the exact 
squared matrix elements, coherently including also interferences among diagrams. The hnal 
step to obtain |A4^‘^p of Eq. (2.4) from Eq. (2.8) is to replace the photon energy spectrum 
with the Altarelli-Parisi distribution for a better treatment of hard collinear radiation. 

®For the sake of simplicity, we consider the decay H —>■ 2e2p + np, the generalization to 4e or 4/r being 
straightforward. 
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Equation (2.4), with the building blocks described above, can then finally be improved 
to include exact NLO corrections according to our master formula: 

1 oo / n \ 

drr'h”* = ^Fsvn({p),e)^_ l'U™({p}.{U)l"d<5„({p},{t)) 

n=0 \z=l / 


Fsv = 1 + 


- dr^,5’" 


sv 


dr^o 


Fha = 1 + 






(2.9) 


The correction factors Fsv and Fn^i carry the information of the exact NLO calcu¬ 
lation: dTgy*^ is the sum of the virtual corrections of Eq. (2.2) and Eq. (2.3) and soft 
real correction given by the first line of Eq. (2.1), dr^y’“ is its PS approximation, i.e. the 
0{a) term without any real hard photon of Eq. (2.4), is the exact one-photon 

bremsstrahhmg amplitude and A4{^‘^’“(/cj) is its PS approximation. 

We want to remark that Fsv and Fj{,i are by construction free of collinear and/or 
infrared logarithms and that the 0{a) expansion of Eq. (2.9) exactly coincides with the NLO 
calculation, without any double counting. Furthermore, Eq. (2.9) is still fully differential 
in the final state momenta and can be conveniently implemented in a MC event generator. 

Finally, we remark that the NLO virtual and real corrections used in Fsv and Fn^i are 
strictly defined only for 0 or 1 photon, while in Eq. (2.9) they are used also when there 
are additional photons: this requires a mapping of the n-photons phase space to 0 or 1 
photon phase space. The mapping is implemented in close analogy to the one described in 
appendix A.2 of Ref. [29], and here we do not discuss it in further detail. 


3 Numerical results 


In the present Section, we show and discuss the numerical results provided by our calcula¬ 
tion, as obtained with the new tool Hto41. First, we show some tuned comparisons with 
the predictions of the reference code Prophecy4f at the level of NLO electroweak correc¬ 
tions. Then, we present our best predictions for various observables at NLOPS electroweak 
accuracy, as well as for Higgs production and decay in the presence of NLO QCD and 
electroweak corrections matched to PS. 

The results presented here are obtained using Prophecy4f v2.0.^ In both codes, we 
use the following set of input parameters 


a(0) = 1/137.03599911 
Mw = 80.398 GeV 
me = 510.99892 KeV 
mu = 190 MeV 
md = 190 MeV 


= 1.16637 lO-*^ GeV-2 
= 2.141 GeV 
mf, = 105.658369 MeV 
me = 1.4 GeV 
ms = 190 MeV 


Mz = 91.1876 GeV 
Tz = 2.4952 GeV 
mr = 1.77684 GeV 
"mtop = 172.5 GeV 
mb = 4.75 GeV 


Table 1. Values of the input parameters used in the numerical calculations. 


^Available at http: //omnibus .uni-freiburg.de/~sd565/programs/prophecy4f/prophecy4f .html 
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The Mz,w and Tz,w are the running-width PDG values which have to be converted 
to the fixed-width scheme adopted here through, for example, the relations of Eq. (7.2) of 
Ref. [35]. As we work in the scheme, for the electromagnetic coupling constant we use 
the expression 

\/2 Gf, sin^ Ow ^ 

ao, = -- (3.1) 

with sin^ 6\y = 1 — /Aff, in the calculation of the LO width and NLO weak corrections, 
while we use q;(0) for the coupling of the photon to the external charged particles.^ The 
top-quark width is set to the LO prediction in the SM, and a fixed width is employed in all 
the resonant propagators in the framework of the complex mass scheme. 

3.1 NLO electroweak corrections: comparisons to Prophecy4f 

A sample of the Prophecy4f vs. Hto41 comparisons at NLO electroweak accuracy is shown 
in Tab. 2 and in Figs. 1-3, in order to check the technical accuracy of our predictions in its 
different aspects sketched in Sect. 2.1. Generally speaking, we observe very good agreement 
between our predictions and the independent results of Prophecy4f. 

In Fig. 1 we show the comparison for the NLO width in the leptonic decay channels 
H —)• 2e2p and H —>■ 4p,® as a function of the Higgs mass in the range [125, 400] GeV, 
together with the relative contribution due to the NLO electroweak corrections where the 
effect of mass thresholds present in the loop computation is particularly visible. As can 
be seen, the two calculations perfectly agree. For the sake of clarity and completeness, we 
quote in Tab. 2 the predictions of the two codes for the decay channels H —?■ 2e2p and 
iL —)• 4// for three specific values of the Higgs mass: the level of agreement is within the 
statistical numerical uncertainty which is well below the 0.1% accuracy. 


Mjy/Final State 

Prophecy4f 

Hto41 

125 GeV/iL ^ 2e2// 

0.24151(8) 

0.24165(2) 

140 GeY/H 2e2// 

1.2672(2) 

1.2667(1) 

200 GeY/H 2e2// 

825.9(1) 

825.8(1) 

125 GeY/H 4// 

0.13324(2) 

0.13325(2) 

140 GeY/H 4/x 

0.6713(1) 

0.6711(1) 

200 GeY/H 4/x 

413.02(7) 

412.98(2) 


Table 2. Comparison between the NLO electroweak predictions of Prophecy4f and our calculation 
(Hto4l) for the Higgs decay width (in KeV), for different values of the Higgs mass and final states. 
The numbers in parenthesis are the statistical uncertainty on the last digit due to MC integration. 

In Fig. 2 a comparison between Prophecy4f and Hto41 is shown for the e"’“e“ invariant 
mass (in the Higgs rest frame), in the range [60, 100] GeV (upper plot) and in the range 

®This value is used for all the numerical results shown in the following, with the exception of the 
comparisons with Prophecy4f, where we use aa^ everywhere, to be consistent with the default choice of 
Prophecy4f. 

^Analogous results are valid in the iif —>■ 4e channel, which coincides for the integrated partial width 
with the 4/r final state (apart from negligible mass effects). 
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Figure 1. Comparison between the NLO electroweak calculation of Prophecy4f and Hto41 for 
the decay widths H —>■ 2e2/i (upper plot) and H ^ Afj. (lower plot), as a function of the Higgs mass 
in the range [125, 400] GeV. For each plot, upper panel: absolute predictions in KeV; lower panel: 
ratio between LO width and NLO EW corrected width. 

[85, 95] GeV (lower plot). The results refer to the decay channel H —>■ 2e2^ for Mh = 
125 GeV. Also in this case, the agreement between the two codes is remarkable, in spite of 
the large effect due to the radiative corrections'. Actually, at and above the peak of the 
electron-pair invariant mass distribution the corrections are of the order of 30%, while for 
Mg+e- below Mz they can reach 50%. The lowering of the peak and the raising of a tail can 
be mainly ascribed to the photon radiation off the leptons, as typical final-state radiation 

^For simplicity, in the present Section we provide results for bare electrons only, i.e. in the absence of 
lepton-photon recombination effects. 


- 9 - 
















Me+e- (GeV) 


Figure 2. Comparison between the NLO electroweak calculation of Prophecy4f and Hto41 for 
the e+e” invariant mass (in the Higgs rest frame), in the range [60, 100] GeV (upper plot) and in 
the range [85, 95[ GeV (lower plot). Predictions for the decay H —>■ 2e2^ at Mh = 125 GeV. Upper 
panels: absolute predictions for dV /lower panels: relative effect of the NLO corrections. 


(FSR) effect observed around the peak of resonant processes [33, 34, 55, 56]. 

A further comparison is given in Fig. 3 for the distribution of the angle between the 
decay planes of the virtual Z bosons in the H rest frame for the channels H —?■ 2e2// (upper 
plot) and H ^ Afj, (lower plot) for Mh = 125 GeV, which is the observable of main interest 
for spin-parity assignment. For the (p angle we use the definition 


(ki2 X ki) • (ki2 X ks) 

^ = “Ti-TTTi- 

|ki2 X ki||ki2 X k3| 

sgn(sin(/)) = sgn{ki 2 • [(ki 2 x ki) x (ki 2 x ka)]} 


(3.2) 

(3.3) 
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Figure 3. Comparison between the NLO electroweak calculation of Prophecy4f and Hto41 for 
the (j> angle distribution (in the Higgs rest frame) for the decay channels H —>■ 2e2/i (upper plot) 
and H —>■ Afj, (lower plot) at Mh = 125 GeV. Upper panels: absolute predictions in GeV/deg; lower 
panels: relative effect of the NLO corrections. 


where ki 2 = ki + k 2 and ki, k 2 , ks, k 4 are the three-momenta of the final-state leptons. 

Again the predictions of the two codes nicely agree. The contribution of the NLO 
corrections is particularly visible at the edges of the distribution, where it can reach the 5% 
level for both the decay channels. 

3.2 Predictions at NLOPS electroweak accuracy 

Some illustrative results obtained according to a number of variants of the theoretical 
approach described in Sect. 2.2 are given in Figs. 4-6. In order to disentangle the impact 
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of the different sources of correction, we consider the results obtained according to the 
following levels of accuracy: 

1. the pure PS approximation for the decay width as in Eq. (2.4), associated to multiple 
photon emission in the soft/collinear limit; 

2. the 0{a) truncated approximation of Eq. (2.4), describing one photon radiation in 
the PS framework; 

3. the complete NLO electroweak calculation; 

4. the NLO QED calculation, given by the gauge-invariant subset of electromagnetic 
contributions within the full set of electroweak corrections; 

5. the NLO electroweak corrections matched to the QED PS, as in Eq. (2.9); 

6 . the NLO QED corrections matched to the QED PS, i.e. the QED gauge-invariant 
realization of Eq. (2.9). 

The comparison between approximations 1. and 2. is useful to quantify the higher- 
orders contribution due to photon emission beyond 0(a), while the difference between 
options 3. and 4. is a measure of pure weak loop corrections, the difference between 
approximations 2. and 3. is an estimate of non-logarithmic 0{a) QED terms plus pure weak 
loop corrections. The comparison between approximations 3. and 5., as well as between 4. 
and 6., allows us to check that the NLOPS matching procedure correctly preserves the effect 
of QED exponentiation as given by the difference between options 1. and 2. Moreover, the 
results of 1. vs. those of 5. and of 3. vs. those of 5. provide an estimate of the accuracy 
of the predictions available in the literature for Higgs physics at the LHC, in particular of 
of the process-independent, widely used code PHOTOS [57], which describes multiple photon 
emission but does not include exact NLO electroweak corrections, and of Prophecy4f, that 
does not take into account the contribution of exclusive QED exponentiation. 

In Fig. 4 we show the relative contribution of the different theoretical approximations 
discussed above for the e'’'e“ (upper plot) and ii~ (lower plot) invariant mass in the 
Higgs rest frame, in the range [85, 95] GeV. The results refer to the process H —)> 2e2/r 
for Mfj = 125 GeV, according to a bare lepton definition. By inspection of Fig. 4 we can 
draw the following conclusions: the NLO corrections to the lepton invariant masses are 
quite large, since they amount to about 50% (30%) to the e"*“e“ (/x^/i”) invariant mass 
below the peak and about 30% (20%) at and above it. They are largely dominated by 
the enhanced leading logarithmic contributions of QED nature oc alog(M^/m|), as can be 
inferred from the comparison between the results of the pure 0{a) PS algorithm and those 
of the NLO QED/electroweak calculations. From this comparison, one can also conclude 
that the 0{a) non-logarithmic QED terms contribute at the some per cent level, both for 
the e'’'e“ and invariant mass, whereas the pure weak loops have a much smaller 

effect, not exceeding the 1% level. 

The large impact of NLO QED corrections, which significantly modify the shape of 
the invariant mass distribution, translates in a relevant contribution due to higher-order 
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Figure 4. Relative contribution of the QED/electroweak corrections to the e“''e“ (upper plot) and 
(lower plot) invariant mass in the Higgs rest frame, in the range [85, 95] GeV. Predictions 
for the decay H —)■ 2e2/x at Mh = 125 GeV. The theoretical approximations corresponding to the 
different lines are explained in the text. 


photonic corrections. Mnltiple photon emission is of the order of 10% for the e'’'e“ hnal- 
state and at the level of some per cents for the case, as a conseqnence of the different 

magnitnde of the lepton-photon collinear logarithm. It can also be noticed that QED 
exponentiation rednces the impact of NLO corrections and that the NLOPS matching 
correctly preserves the size of mnltiple photon emission. 

Qnite different conclnsions derive from the analysis of Fig. 5, which shows the relative 
corrections of the different theoretical recipes on the (p angle distribntion for the H —?• 2e2^ 
and —7- 4^ decays. For snch an observable, the pnre 0(a) PS approximation signihcantly 
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Figure 5. The same as Fig. 4 for the (j) angle distribution in the decay channels H —)• 2e2/i (upper 
plot) and H ^ An (lower plot) at Mh = 125 GeV. 

underestimates the contribution of NLO EW corrections for (j) close to 0° and 360°, while 
it provides an overestimate around 180°. Actually, it can be noticed that the 4> angle 
distribution receives a non-negligible contribution from fixed-order non-logarithmic terms 
and that, more importantly, is particularly sensitive to pure weak corrections, which set the 
correct overall size and shape of the radiative corrections. On the other hand, the effect of 
QED exponentiation is moderate, varying between a few per mille to about 1%. 

For completeness, we show in Fig. 6 results for the invariant mass of the e"’'e“ pair 
and the (j) angle distribution (for the process H —)• 2e2^) under the more realistic exper¬ 
imental condition of calorimetric or recombined electrons and positrons. In this case, we 
replace the three-momentum of the with the effective momentum p = Pg± -|- p^ for each 
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Figure 6. Relative contribution of the QED/electroweak corrections to the e+e” invariant mass 
(upper plot) and the (j) angle distribution (lower plot) for recombined electrons and positrons. 
Predictions for the decay H —)■ 2e2/r at Mh = 125 GeV in the Higgs rest frame. The theoretical 
approximations corresponding to the different lines are explained in the text. 

photon satisfying the condition < 0.1, as typically done 

by LHC experiments, where A(f)f,±^ is the lepton-photon separation angle in the transverse 
plane. As can be seen from Fig. 6 in comparison to Fig. 4 and Fig. 5, the contribution 
of the radiative corrections is largely reduced, as expected, when switching from bare to 
recombined electrons/positrons. For the e'’'e“ invariant mass, the corrections are reduced 
by about a factor of three, almost independently of the considered theoretical approxima¬ 
tion, and preserve their shape. However, non-negligible corrections still remain under the 
calorimetric condition, of about +15% in the left tail of the invariant mass and of the order 
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of -10% around the peak of the distribution, when considering the most accurate matched 
predictions. In comparison to the case of bare electrons, the effect of QED exponentiation 
for dressed electrons reduces to about 1% in the tail and at the per mille level at and above 
the peak. 

More interestingly, the QED relative corrections to the (p angle distribution are sub¬ 
stantially modified both in size and shape by the recombination effects, whereas the full 
electroweak predictions receive a slight size reduction and a less pronounced shape modih- 
cation. In particular, we checked through detailed numerical inspections that the especially 
visible difference in shape between the pure PS and the diagrammatic QED predictions is 
of virtual origin and has to be ascribed to the QED pentagons, which are exactly included 
in the Feynman diagram calculation and only (crudely) approximated in the soft/collinear 
limit in the PS calculation. To some extent, we expect that the rich angular correlations 
introduced by pentagon diagrams is only poorly reproduced by the PS approximation. All 
in all, the results shown in the lower plot of Fig. 6 reinforce the already noticed particularly 
relevant role played by loop contributions with complex topology, both of QED and weak 
nature, to obtain reliable predictions for the (p angle observable. 

To summarize, the main conclusion of this Section is that both NLO electroweak and 
higher-order QED corrections, as well as their combination, are relevant for reliable simula¬ 
tions of the most important observables considered in precision studies of the Higgs sector 
at the LHC. 

3.3 Interface to POWHEG: results for production and decay 

In order to facilitate phenomenological studies of Higgs boson production and decay in the 
presence of both QCD and electroweak contributions, we have implemented an interface 
which allows to use our code in association with any event generator describing Higgs pro¬ 
duction. In Figs. 7-9 we show a sample of illustrative results obtained by interfacing Hto41 
with POWHEG [42] for the simulation of Higgs boson production in gluon-gluon fusion. We 
use the POWHEG version with NLOPS accuracy in QCD [58] from the POWHEG BOX frame¬ 
work [59] and we consider Higgs production in proton-proton collisions at a c.m. energy of 
8 TeV®. The events generated by POWHEG are interfaced to Hto41 according to the following 
procedure: 

• generate unweighted events for the process pp —)• H[+j) in the Les Houches format 
using POWHEG, where H is an on-shell Higgs boson and j stands for the extra parton 
of the NLO QCD calculation; 

• the Les Houches hie is read event by event by Hto41 and the particles momenta are 
stored in the generic common block structure introduced in Ref. [60]; 

• each event is decayed into the selected channel in the H rest frame, using Hto41. 
After boosting the decay products back to the laboratory frame, the events including 
production and decay are written in a hie in the Les Houches format. 

® However, as we are interested to study the relative impact of electroweak corrections dominated by 
contributions of the kind a" log"'(M|/m|), the results shown in the following are in practice independent 
of the c.m. energy. 
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Figure 7. Comparison between the results obtained with POWHEG + Hto41 (Born) + PYTHIA v6 
(red dashed line) and POWHEG + Hto41 (NLOPS) + PYTHIA v6 (blue solid and black dash-dotted 
lines) for the transverse momentum (upper plot) and rapidity (lower plot) of the Higgs boson. In 
the lower panels the relative contribution of NLOPS electroweak corrections for bare (solid line) 
and recombined (dash-dotted line) leptons is shown. 


The Les Houches file can be finally passed to a shower event generator for QCD shower¬ 
ing and hadronization. In our examples we use PYTHIA v6.4 [61] as QCD PS. According to 
the above procedure, the pp ^ H ^ M process is treated in narrow width approximation, 
as it is the case for a 125 GeV Higgs boson, and factorized in on-shell Higgs production and 
decay. 

In our analysis we consider, for definiteness, the decay channel H —)• 2e2/i and the 
following observables: the transverse momentum pi^ and rapidity yn of the Higgs boson 


17 - 





































cos 6* 


Figure 8. The same as in Fig. 7 for the invariant mass of the subleading lepton pair (upper plot) 
and the cosine of the angle of the leading lepton pair in the four-lepton rest frame with respect to 
the beam axis (lower plot). 


(Fig. 7), the invariant mass of the subleading lepton pairs and the magnitude of the cosine 
of the decay angle of the leading lepton pair in the four-lepton rest frame with respect 
to the beam axis |cos0*| (Fig. 8). The leading pair is defined as the lepton pair with 
invariant mass closest to the Z boson mass and its angle is obtained by summing the 
three-momenta of the two leptons. For the POWHEG calculation of Higgs production in gluon 
fusion, we use the PDF set MSTW2008nlo68cl [62] with factorization/renormalization scale 
IJ'R = fJ-F = Mh- The values of the other input parameters are the same as the ones given 
in Tab. 1. The results shown in the following refer to a sample of 1, 2 • 10® unweighted events 
and to the same selection cuts adopted in Ref. [23] and correspond to bare (solid line) and 
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Figure 9. Distribution of the transverse momentum of the hardest photon (upper plot) and the 
angular separation between the hardest photon and the closest lepton (lower plot) obtained using 
POWHEG + Hto41 (NLOPS) + PYTHIA v6. The minimum photon energy is if™™ = 1 GeV. 

recombined (black dash-dotted line) leptons. In the latter case, we recombine photons with 
both electrons and muons, in analogy to the selection criteria adopted in the experimental 
study of Ref. [23], if the condition < 0.1 is satisfied. 

In Fig. 7 and Fig. 8 we show the comparison between the predictions obtained using 
POWHEG interfaced to our code at LO and NLOPS electroweak accuracy. It can be noticed 
that the contribution due to NLOPS electroweak corrections is almost flat and of about 
— 15 (—5)% for , uh and | cos0*| when considering bare (recombined) leptons, while the 
invariant mass of the subleading lepton pairs receives a varying correction of size between 
—20 (—10)% and —10 (—5)% for bare (calorimetric) leptons, respectively. 
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In Fig. 9 we show the results for two observables which are fully exclusive over QED 
radiation and which can be easily treated in our approach. The results correspond to the 
process pp —?■ —?• 2e2/i + nj, with FI™™ = 1 GeV, for which we show the transverse 

momentum of the hardest photon and the angular separation between the hardest photon 
and the closest lepton, that exhibit the expected features of photon emission in radiative 
events. 


4 Conclusions 

In this work we have presented a precision calculation of the SM Higgs boson decay into 
four charged leptons, in view of improved measurements of the properties of the Higgs 
particle at the LHC Run H. Our approach is based on the computation of the full one-loop 
electroweak corrections supplemented with the contribution of multiple photon emission 
taken into account according to a fully exclusive QED PS algorithm. Our results, which 
have a NLOPS electroweak accuracy, are available in the form of a new event generator, 
Hto41, that can be easily interfaced to any QCD program simulating Higgs production. 

We have cross-checked our NLO electroweak corrections against the predictions of the 
reference code Prophecy4f and found perfect agreement. We have also shown that both 
NLO electroweak and higher-order QED corrections, as well as their interplay, are necessary 
for actually precise simulations of the variety of observables involved in Higgs physics at 
the LHC. This provides the main novel theoretical feature of our work, which goes beyond 
the presently available results limited to the fixed-order approximation or to a leading loga¬ 
rithmic QED modeling. The second relevant aspect is given by the possibility of interfacing 
Hto41 to any generator describing Higgs boson production, thus allowing simulations of 
Higgs production and decay in the presence of higher-order QCD and electroweak correc¬ 
tions. In this respect, we have shown some illustrative results obtained in terms of the 
combined usage of POWHEG and Hto41. 

Our results can find application in precision measurements of the Higgs boson mass, 
spin-parity determination and tests of the SM at the level of differential cross sections in 
the future run of the LHC. They can be generalized to other processes yielding four leptons 
in hadronic collisions, like e.g. pp ^ H ^ —>■ 2i2v or pp —)• —)• Al. 
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A Phase space parameterisation and integration 


The 4 + n bodies phase space as in Eq. (2.5) is integrated according to standard multi¬ 
channel MC techniques, combined with importance sampling to reduce the variance of the 
integral and help event generation®. The first step is to generate a photon multiplicity n 
and associate ni ( 712 ) photons to the electron (muon) current {n\ + n 2 = n), dehning the 
channel of the multi-channel integration. The phase space is then conveniently split into 
two decaying objects to follow the Z propagators, namely 

d4>(Pff;pi,--- - ,kn) = {2TifdQ%^dQ\^<l^{PH\Pzi,Pz2)^ 

d4>(Pz,;pi,p2,^i,--- ,^ni) d4>(Pz2;p3,p4,A:ni+i,--- ,kni+n2) (AT) 

where P^. (P|, = are the momenta of the virtual Z bosons. 

We refrain from writing explicitly the simple 1 —>■ 2 decay phase spaces of Eq. (A.l) 
and we focus instead on the case where at least one photon is present. As discussed in 
appendix A.3 of Ref. [29], an efficient sampling of photons collinear to hnal state leptons 
is a non trivial task, because the directions of the leptons are known only after all the 
momenta are generated. In Ref. [29] we adopted a solution based on a properly chosen 
multi-channel strategy. Here we adopt a different and elegant solution, which consists in 
writing the phase space in the frame where the leptons are back-to-back, i.e. pa = —pb (see 
for example [63-65]). 

Omitting overall numerical factors for brevity, the building block we are interested in 
is 


d$(P;Pa,Pfe,fei,--- ,kr) = (p-Pa-Pb-'^k}\ 

V i=l J Pb O 

= 6^^\p - Q - K) 6<l> 

where we dehned Q = Pa + Pb, K = contains the inhnitesimal phase space 

element divided by the hnal state particle energies. It is usually understood that all the 
variables are expressed in the frame where P is at rest, but we want to express them where 
Q is at rest. In order to do that, the previous equation can be further manipulated by 
inserting the following identities 

d'^Q (Q -pa-pb) = 1 

ds'6{Q^ - s') = 1 

d^P (P) <5(p0 - ^^) = 2\/i d^P (P) <5(p2 - s) = 1 (A.2) 

which help to make explicit the Lorentz invariance of the phase space element. 

®Here we consider only the decay H —> 2e2fi, the generalization to 4 identical leptons being straightfor¬ 
ward. 
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With the help of Eq. (A.2) and appropriately rearranging the terms, we can write 


- = <5$ d^Q6^^\Q - Pa - Pb)S^^HP) 

2^s d^P5(^)(P -Q- K)6{P^ - s)ds'dQ°(5(Q2 _ g') = 

= 6<!>^ -Pa- Pb)S{{Q + Kf - s)d3Qds' = 

. - P.-vmQ + Kf - n 


Pa Pb 


k? 


-•- “T I—. 2=1 * 


(A.3) 


In the cascade of identities (A.3) we used the result (PQ = ( 3 ^ 5)2 (see [64]) and we 

made use of Lorentz invariance. In the last identity it is understood that all the variables 
are expressed in the frame where Q = pa + Pb is at rest and s' = s = /3a is the 

speed of particle a and dfla = dcos0ad(/>a- The big advantage of the last equation is that 
the lepton momenta Pa and pf, lie on the same direction defined by cos 9a and (/a, hence all 
photons can be generated along this direction to sample the collinear singularities. Once all 
particle momenta are generated, they can be boosted back to the rest frame of the decaying 
Higgs boson. 

One last remark concerns the integration limits of the phase space. As mentioned in 
Sect. 2.2, photon energies should be generated larger than the infrared cut-off e in the Higgs 
frame, which is a non Lorentz invariant cut. Since the minimum photon energy can not be 
determined a priori in the frame where Q is at rest (because Q itself depends on the photons 
momenta), we decide to generate photon energies starting from 0 to cover the whole phase 
space and then, once boosted back, cut the event if a photon enegy falls below e. Finally, in 
order to flatten the infrared divergence, we choose to sample the photon energies according 
to the function 


f{ui) oc 



u}>e' 
oj <e' 


where e' is a guessed (and tuned for efficiency) 


minimum energy. 


References 

[1] ATLAS Collaboration, G. Aad et al., Observation of a new particle in the search for the 
Standard Model Higgs boson with the ATLAS detector at the LHC, Phys.Lett. B716 (2012) 
1-29, [arXiv; 1207.7214]. 

[2] CMS Collaboration, S. Chatrchyan et al., Observation of a new boson at a mass of 125 GeV 
with the CMS experiment at the LHC, Phys.Lett. B716 (2012) 30-61, [arXiv; 1207.7235]. 

[3] F. Englert and R. Brout, Broken symmetry and the mass of gauge vector mesons, Phys. Rev. 
Lett. 13 (Aug, 1964) 321-323. 

[4] P. Higgs, Broken symmetries, massless particles and gauge fields. Physics Letters 12 (1964), 
no. 2 132 - 133. 


22 - 







[5] P. W. Higgs, Broken symmetries and the masses of gauge bosons, Phys. Rev. Lett. 13 (Oct, 
1964) 508-509. 

[6] G. S. Guralnik, G. R. Hagen, and T. W. B. Kibble, Global conservation laws and massless 
particles, Phys. Rev. Lett. 13 (Nov, 1964) 585-587. 

[7] P. W. Higgs, Spontaneous symmetry breakdown without massless bosons, Phys. Rev. 145 
(May, 1966) 1156-1163. 

[8] T. W. B. Kibble, Symmetry breaking in non-abelian gauge theories, Phys. Rev. 155 (Mar, 
1967) 1554-1561. 

[9] ATLAS Gollaboration, G. Aad et al.. Measurements of Higgs boson production and 
couplings in the four-lepton channel in pp collisions at center-of-mass energies of 7 and 8 
TeV with the ATLAS detector, Phys.Rev. D91 (2015), no. 1 012006, [arXiv; 1408.5191]. 

[10] ATLAS Gollaboration, G. Aad et al.. Measurement of Higgs boson production in the 
diphoton decay channel in pp collisions at center-of-mass energies of 7 and 8 TeV with the 
ATLAS detector, Phys.Rev. D90 (2014), no. 11 112015, [arXiv; 1408.7084]. 

[11] ATLAS Gollaboration, G. Aad et al.. Measurement of the Higgs boson mass from the 

H —> 77 and H —>■ ZZ* —>■ 4£ channels with the ATLAS detector using 25 fb~^ of pp collision 
data, Phys.Rev. D90 (2014), no. 5 052004, [eirXiv: 1406.3827]. 

[12] CMS Gollaboration, V. Khachatryan et al.. Precise determination of the mass of the Higgs 
boson and tests of compatibility of its couplings with the standard model predictions using 
proton collisions at 7 and 8 TeV, arXiv: 1412.8662. 

[13] CMS Collaboration, V. Khachatryan et al.. Observation of the diphoton decay of the Higgs 
boson and measurement of its properties, Eur.Phys.J. C74 (2014), no. 10 3076, 

[arXiv: 1407.0558]. 

[14] CMS Collaboration, S. Chatrchyan et al.. Measurement of the properties of a Higgs boson in 
the four-lepton final state, Phys.Rev. D89 (2014) 092007, [arXiv: 1312.5353]. 

[15] ATLAS Collaboration, G. Aad et al.. Measurements of Higgs boson production and 
couplings in diboson final states with the ATLAS detector at the LHC, Phys.Lett. B726 
(2013) 88-119, [arXiv: 1307.1427]. 

[16] ATLAS Collaboration, G. Aad et al.. Evidence for the Higgs-boson Yukawa coupling to tau 
leptons with the ATLAS detector, arXiv: 1501.0494. 

[17] CMS Collaboration, S. Chatrchyan et al.. Evidence for the direct decay of the 125 GeV 
Higgs boson to fermions. Nature Phys. 10 (2014) [arXiv: 1401.6527]. 

[18] CMS Collaboration, V. Khachatryan et al.. Constraints on the spin-parity and anomalous 
HVV couplings of the Higgs boson in proton collisions at 7 and 8 TeV, arXiv: 1411.3441. 

[19] ATLAS Collaboration, G. Aad et al.. Evidence for the spin-0 nature of the Higgs boson 
using ATLAS data, Phys.Lett. B726 (2013) 120-144, [arXiv: 1307.1432]. 

[20] LHC Higgs Cross Section Working Group Collaboration, S. Dittmaier et al.. Handbook 
of LHC Higgs Cross Sections: 1. Inclusive Observables, eirXiv: 1101.0593. 

[21] S. Dittmaier, C. Mariotti, G. Passarino, R. Tanaka, et al.. Handbook of LHC Higgs Cross 
Sections: 2. Differential Distributions, arXiv: 1201.3084. 

[22] LHC Higgs Cross Section Working Group Collaboration, S. Heinemeyer et al.. 
Handbook of LHC Higgs Cross Sections: 3. Higgs Properties, arXiv: 1307.1347. 


23 - 


[23] ATLAS Collaboration, G. Aad et al., Fiducial and differential cross sections of Higgs boson 
production measured in the four-lepton decay channel in pp collisions at y/s=8 TeV with the 
ATLAS detector, Phys.Lett. B738 (2014) 234-253, [arXiv: 1408.3226]. 

[24] ATLAS Collaboration, G. Aad et al.. Determination of the off-shell Higgs boson signal 
strength in the high-mass ZZ and WW final states with the ATLAS detector, 

arXiv;1503.0106. 

[25] CMS Collaboration, V. Khachatryan et al.. Constraints on the Higgs boson width from 
off-shell production and decay to Z-boson pairs, Phys.Lett. B736 (2014) 64, 

[arXiv: 1405.3455]. 

[26] C. Carloni Calame, C. Lunardini, G. Montagna, O. Nicrosini, and F. Piccinini, Large angle 
Bhabha scattering and luminosity at flavor factories, Nucl.Phys. B584 (2000) 459-479, 
]hep-ph/0003268] . 

[27] C. Carloni Calame, G. Montagna, O. Nicrosini, and M. Treccani, Higher order QED 
corrections to W boson mass determination at hadron colliders, Phys.Rev. D69 (2004) 
037301, ]hep-ph/0303102]. 

[28] C. M. Carloni Calame, G. Montagna, O. Nicrosini, and M. Treccani, Multiple photon 
corrections to the neutral-current Drell-Yan process, JHEP 0505 (2005) 019, 
]hep-ph/0502218]. 

[29] G. Balossini, C. M. Carloni Calame, G. Montagna, O. Nicrosini, and F. Piccinini, Matching 
perturbative and parton shower corrections to Bhabha process at flavour factories, Nucl.Phys. 
B758 (2006) 227-253, ]hep-ph/0607181]. 

[30] G. Balossini, C. Bignamini, C. C. Calame, G. Montagna, O. Nicrosini, et al.. Photon pair 
production at flavour factories with per mille aceuracy, Phys.Lett. B663 (2008) 209-213, 

[arXiv :0801. 3360] . 

[31] Working Group on Radiative Corrections and Monte Carlo Generators for Low 
Energies Collaboration, S. Actis et al.. Quest for precision in hadronic cross sections at low 
energy: Monte Carlo tools vs. experimental data, Eur.Phys.J. C66 (2010) 585-686, 

[arXiv: 0912. 0749]. 

[32] C. Carloni Calame, G. Montagna, O. Nicrosini, and F. Piccinini, The BABAYACA event 
generator, Nucl.Phys.Proc.Suppl. 131 (2004) 48-55, ]hep-ph/0312014]. 

[33] C. Carloni Calame, G. Montagna, O. Nicrosini, and A. Vicini, Precision electroweak 
calculation of the charged current Drell-Yan process, JHEP 0612 (2006) 016, 
]hep-ph/0609170]. 

[34] C. Carloni Calame, G. Montagna, O. Nicrosini, and A. Vicini, Precision electroweak 
calculation of the production of a high transverse-momentum lepton pair at hadron colliders, 
JHEP 0710 (2007) 109, [arXiv: 0710. 1722]. 

[35] A. Bredenstein, A. Denner, S. Dittmaier, and M. Weber, Precise predietions for the 

Higgs-boson decay H WWjZZ —>■ 4 leptons, Phys.Rev. D74 (2006) 013004, 

]hep-ph/0604011]. 

[36] A. Bredenstein, A. Denner, S. Dittmaier, and M. Weber, Radiative corrections to the 
semileptonic and hadronic Higgs-boson decays H —>■ WWjZZ —>■ Jf. fermions, JHEP 0702 
(2007) 080, ]hep-ph/0611234]. 


24 - 


[37] A. Bredenstein, A. Denner, S. Dittmaier, and M. Weber, Precision calculations for the Higgs 
decays H —>■ ZZjWW —>■ 4 leptons, Nucl.Phys.Proc.Suppl. 160 (2006) 131-135, 
[hep-ph/0607060]. 

[38] A. Bredenstein, A. Denner, S. Dittmaier, and M. Weber, Precision calculations for 

H WWjZZ —>■ 4 fermions with PROPHECY4^f arXiv;0708.4123. 

[39] F. Piccinini, C. Carloni Calame, G. Montagna, O. Nicrosini, M. Moretti, et al., QED 
corrections to Higgs boson decay into four leptons at the LHC, PoS HEP2005 (2006) 307. 

[40] C. Carloni Calame, M. Moretti, G. Montagna, O. Nicrosini, F. Piccinini, et al.. Impact of 
QED corrections to Higgs decay into four leptons at the LHC, Nucl.Phys.Proc.Suppl. 157 
(2006) 73-77, [hep-ph/0604033[. 

[41] P. Nason, A New method for combining NLO QCD with shower Monte Carlo algorithms, 
JHEP 0411 (2004) 040, [hep-ph/0409146[. 

[42] S. Frixione, P. Nason, and C. Oleari, Matching NLO QCD computations with Parton Shower 
simulations: the POWHEC method, JHEP 0711 (2007) 070, [corXiv;0709.2092]. 

[43] G. ’t Hooft and M. J. Veltman, Scalar One Loop Integrals, Nucl.Phys. B153 (1979) 365-401. 

[44] G. Montagna, O. Nicrosini, and G. Passarino, Analytic final state corrections to e^e“ — >■ // 
with realistic cuts, Phys.Lett. B309 (1993) 436-442. 

[45] A. Denner, S. Dittmaier, M. Roth, and L. Wieders, Electroweak corrections to 
charged-current e+e“ —>■ fermion processes: Technical details and further results, 

Nucl.Phys. B724 (2005) 247-294, [hep-ph/0505042[. 

[46] A. Denner and S. Dittmaier, The Complex-mass scheme for perturbative calculations with 
unstable particles, Nucl.Phys.Proc.Suppl. 160 (2006) 22-26, [hep-ph/0605312[. 

[47] J. Vermaseren, New features of FORM, math-ph/0010025. 

[48] T. Hahn and M. Perez-Victoria, Automatized one loop calculations in four-dimensions and 
D-dimensions, Comput.Phys.Commun. 118 (1999) 153-165, [hep-ph/9807565[. 

[49] T. Hahn, Automatic loop calculations with FeynArts, FormCalc, and LoopTools, 
Nucl.Phys.Proc.Suppl. 89 (2000) 231-236, [hep-ph/0005029[. 

[50] A. Denner and S. Dittmaier, Reduction of one loop tensor five point integrals, Nucl.Phys. 
B658 (2003) 175-202, [hep-ph/0212259[. 

[51] A. Denner and S. Dittmaier, Scalar one-loop 4-point integrals, Nucl.Phys. B844 (2011) 
199-242, [arXiv;1005.2076]. 

[52] G. Passarino and M. Veltman, One Loop Corrections for e'^e~ Annihilation Into in 

the Weinberg Model, Nucl.Phys. B160 (1979) 151. 

[53] A. Denner, Techniques for calculation of electroweak radiative corrections at the one loop level 
and results for W physics at LEP-200, Fortsch.Phys. 41 (1993) 307-420, [airXiv: 0709.1075]. 

[54] B. A. Kniehl, Higgs phenomenology at one loop in the standard model, Phys.Rept. 240 (1994) 
211-300. 

[55] L. Barze, G. Montagna, P. Nason, O. Nicrosini, and F. Piccinini, Implementation of 
electroweak corrections in the POWHEC BOX: single W production, JHEP 1204 (2012) 037, 
[arXiv: 1202.0465]. 


- 25 


[56] L. Barze, G. Montagna, P. Nason, O. Nicrosini, F. Piccinini, et al., Neutral current Drell-Yan 
with combined QCD and electroweak corrections in the POWHEG BOX, Eur.Phys.J. C73 
(2013), no. 6 2474, [arXiv; 1302.4606], 

]57] E. Barberio and Z. Was, PHOTOS: A Universal Monte Carlo for QED radiative corrections. 
Version 2.0, Comput.Phys.Commun. 79 (1994) 291-308. 

]58] S. Alioli, P. Nason, C. Oleari, and E. Re, NLO Higgs boson production via gluon fusion 
matched with shower in POWHEG, JHEP 0904 (2009) 002, [cirXiv;0812.0578]. 

]59] S. Alioli, P. Nason, C. Oleari, and E. Re, A general framework for implementing NLO 

calculations in shower Monte Carlo programs: the POWHEG BOX, JHEP 1006 (2010) 043, 
[arXiv; 1002.2581]. 

[60] E. Boos, M. Dobbs, W. Giele, I. Hinchliffe, J. Huston, et al.. Generic user process interface 
for event generators, hep-ph/0109068. 

[61] T. Sjostrand, S. Mrenna, and P. Z. Skands, PYTHIA 6.4 Physics and Manual, JHEP 0605 
(2006) 026, ]hep-ph/0603175]. 

[62] A. Martin, W. Stirling, R. Thorne, and G. Watt, Parton distributions for the LHC, 
Eur.Phys.J. C63 (2009) 189-285, [arXiv:0901.0002]. 

[63] M. Schonherr and F. Krauss, Soft Photon Radiation in Particle Decays in SHERPA, JHEP 
0812 (2008) 018, [arXiv:0810.5071]. 

[64] S. Jadach, Monte-Carlo Methods for High-Energy Phyiscs, Lectures at the Torino School of 
Physics (2001) [http://jadach.home.cern.ch/jadach/]. 

[65] S. Jadach, B. Ward, and Z. Was, The Precision Monte Carlo event generator KK for two 
fermion final states in e^e“ collisions, Comput.Phys.Commun. 130 (2000) 260-325, 
]hep-ph/9912214]. 


26 


